Re-study on the contribution of scalar potential and spectra of cc, bb and bc(bc) families 
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We indicated in our previous work that for QED the role of the scalar potential which appears at the loop level 
is much smaller than that of the vector potential and in fact negligible. But the situation is different for QCD, 
one reason is that the loop effects are more significant because a s is much larger than a, and secondly the 
non-perturbative QCD effects may induce a sizable scalar potential. In this work, we phenomenologically study 
the contribution of the scalar potential to the spectra of charmonia, bottomonia and bc(bc) family. Taking into 
account both vector and scalar potentials, by fitting the well measured charmonia and bottomonia spectra, we 
re-fix the relevant parameters and test them by calculating other states of not only the charmonia, bottomonia, 
but also further the be family. We also consider the Lamb shift of the spectra. 



I. INTRODUCTION 

The potential model has been proposed to evaluate the spec- 
tra of the quantum systems composing of heavy flavors, such 
as charmonia and bottomonia for many years fljj] . The sub- 
ject on heavy quarkonia was thoroughly discussed in an en- 
lightening paper [2]. Even though one calculates the binding 
energy in terms of the non-relativistic Schrodinger equation, 
it is indeed a reasonable framework for such heavy resonant 
states. The main physics, no matter it is induced by the stan- 
dard model (SM) or new physics beyond the SM, is included 
in the potential. In phenomenology, the potential contains 
two pieces, the Coulomb potential which is induced by the 
gluon exchange and the confinement which may come from 
the non-perturbative QCD effects or other sources. It was 
suggested that due to the symmetry consideration [3], for the 
QED case, if the Dirac equation has a higher symmetry degree 
than the corresponding Schrodinger equation, the Coulomb 
potential should have both scalar and vector parts with equal 
fractions. However, for such a combination the hydrogen po- 
tential would not possess the coupling of orbital angular mo- 
mentum and spin (L • S ) and it definitely contradicts to the 
reality. 

In our previous work 0], we analyzed the source of the 
scalar and vector potentials for the QED, namely the vector 
potential is induced by one-photon exchange for the vector- 
like gauge theory QED, while the scalar part must be com- 
ing from the loop effects. Therefore, for the QED case, the 
fraction of the scalar potential is very suppressed because the 
coupling a ~ 1/137 is small and does not have substantial 
contribution to the spectrum of hydrogen. This conclusion 
is consistent with the precise measurement on the spectra of 
hydrogen-like atoms. However, for the QCD case, the situa- 
tion is very different. First because the strong coupling a s is 
much larger than a, the loop suppression is not as much as 
for the QED and sometimes the NLO effect of QCD can even 
exceed the leading order||5l |6[] . Secondly the non-perturbative 
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effects may also change the whole scenario. Thus for QCD, 
the scalar potential may play an important role and its contri- 
bution may especially manifest at values of the spectrum line 
splitting due to the coupling of orbital angular momentum and 
spin. That was also suggested by some aufhors|01. 

Because the non-perturbative QCD effects cannot be reli- 
ably derived from any underlying theory so far, we prefer 
to introduce phenomenological parameters to manifest their 
roles and the parameters are fixed by fitting the well-measured 
data. Indeed the parameters not only contain the contribu- 
tion of the non-perturbative effects, but also that from the 
higher orders of the perturbative effects. In our last work [8], 
we introduced two more parameters to account for the scalar 
piece in the potential and then re-fitted the spectra of char- 
monia. The fitting is indeed improved comparing with that 
including only the vector piece. Thus for further investigat- 
ing the involvement of scalar potential, we use the same strat- 
egy adopted in fH] to evaluate the spectra of the bottomonia 
(bb) and then the various resonant states of the bc(bc) system. 
Namely, we write U(r) — V(r) + S(r) and the vector po- 
tential is V(r) — — c CFCtsjr + dn 2 r, while the scalar one 
is S(r) = —(a — c)CVa s /f + {b — d)n 2 r with a, b, c, d are 
four parameters to be determined. The induced terms, such 
as L • S coupling, Si • S2 coupling and the spin-independent 
corrections etc. are given in the very enlightening work by 
Lucha et al. i^- flTtl . Substituting all the expressions into 
the Schrodinger equation, we may obtain the spectra of heavy 
quark-antiquark systems. 

Besides, the QED theory predicts the Lamb shift which is 
due to the vacuum effects. In QM, it only shifts the S-wave 
spectra because in the non-relativistic limit, it is proportional 
to S(r), but by the quantum field theory, the L / states are 
also affected. In other words, by considering the Lamb shift, 
the positions of the spectra would deviate from that obtained 
without the Lamb shift. For the hadron case, the governing 
theory is QCD which also induces the Lamb shift Il2ll . and 
in this work, we include its contribution. It is noticed that, 
as the phenomenological parameters which are determined by 
fitting data are introduced, all higher order effects should also 
be automatically involved, so it seems that there is no need 
to consider the Lamb shifts which are induced by the loop 
effects. In fact, it is. However, as calculating the form fac- 
tors of the hadronic transition matrix elements or obtaining 
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the parton distribution functions, we always wish to squeeze 
the uncontrollable parts which are not calculable, such as the 
non-perturbative contributions, as small as possible. Simi- 
larly, here we include the NLO or even NNLO corrections 
i.e. the Lamb shifts and re-fit the parameters which are indeed 
not derivable. 

Considering the Lamb shift, there is a byproduct, which 
may be very enlightening for understanding the theory. Obvi- 
ously, only the products ca s and dn in the potential matter, but 
not c, d, a s and k separately. However, for explicitly showing 
the roles of the scalar and vector pieces, we adopt a special 
strategy. 

When the Lamb shift is taken into account, the situation 
may change slightly. In the expression of the Lamb shift, 
there is an ultraviolet divergent term which includes the renor- 
malization scale fi. Meanwhile the running coupling a s also 
depends on the scale. The authors of II 1311 suggest an effec- 
tive method to deal with the divergence and meanwhile fix the 
value of a s (see the text, where we introduce the method in 
some details for the readers' convenience). It is noted that ji 
is a complicated function of a s , quark mass (mj or m c ), and 
the principal quantum number n. Taking a special way to de- 
termine [i which is corresponding to adopting a special renor- 
malization scheme, and considering the dependence of a s on 
fi, one can eventually find the value of a s for a certain flavor 
(6 or c) and a given principal quantum number. For example, 
for T(1S), we have a s = 0.218 for m b = 4.8 GeV and the 
scale-parameter A = 0.2 GeV. Amazingly, this value is quite 
close to that adopted in literature by fitting the T spectra. 

By contrast, the confinement term nr is fully coming from 
the non-perturbative QCD effects and cannot be theoretically 
derived so far. Thus we adopt the value given in lfl5ll . 

With all the inputs, we calculate the spectra of the cc, bb 
and bc(bc) systems. 

This paper is organized as follows. In Sec. [II] and [TTH 
we introduce the generalized Breit-Fermi Hamiltonian and 
the Schrodinger equation for the bb bound states: T(1S), 
Xbo(lP), Xfci(lP), Xm(1P) and T(2S). Then we numerically 
solve the eigen-equations for these bound states and fix the 
parameters as we did for dealing with the charmonia family 
in our previous work. In Sec. |IV| the Lamb shift is taken into 
account and another set of the parameters is given to improve 
our predictions. In Sec.[VJ we give the spectra of the bc(bc) 
mesons. The last section is devoted to our conclusion and dis- 
cussion. 



II. THE GENERALIZED BREIT-FERMI HAMILTONIAN 
AND SCHORDINGER EQUATION 



The generalized Breit-Fermi Hamiltonian was given in 
Refs. JM1 as 



Hi — H sc \ + H si , 

H s d = His + H ss + Hi 
1 
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(le) 



where, V and S stand as the vector and scalar potentials 
and H S [ and H s d represent the spin-independent and spin- 
dependent pieces respectively. For the linear confinement 
piece we adopt the Cornell potential lfl6ll . Thus the total po- 
tential at the lowest order reads 



U(r) = V(r) + S(r) 



-aCp 



bn r 



where, 



V(r) 
Sir)- 



-c CFCt s /r + dn 2 r 

-(a — c)C F a s /r + (b - d)n 2 r 



(2a) 



(2b) 



With the hamiltonian ([T) and the potential ©, one can solve 
the Schrodinger equation 



H^(r) = (H + i?i)*(r) = (E + 2m)#(r). 



(3) 



We can transform the radial wave function into R(x) with 
a dimensionless variable: x — nr, then the reduced radial 
equation is written as 1 

d 2 



dx 2 



u(x) = A(x)u(x), 



(4a) 



where, 



A(x) = -m(E- U(x) - H-i 

1 /- ~ \ 2 
-L{e-U(x) 



1(1 + 1) 



with 

m = m/n, E — E/k, 
Hi = Hi/k, U(x) = U(x)/k. 

In the simplified potential form (O, the approximation 
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P 



m(E-U(r)) 



(4b) 



(4c) 



(4d) 



is used. For the legitimacy of applying this approximation in 
the calculation and the error degree brought up in the numeri- 
cal values are briefly discussed in the appendix. 



H = H Q + Hi + . . . , (la) 

2 

H = ^ + 2m + S(r) + V(r) , (lb) 



1 The standard form of the radial equation can be easily found in Ref. fl7h . 
and the method to make it dimensionless is borrowed from Ref. fT^I . 
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III. THE ENERGY GAP FUNCTION FOR THE bb 
BOTTOMONIA AND THE NUMERICAL RESULTS 
WITHOUT TAKING INTO ACCOUNT THE LAMB SHIFT 

The radial equation (|4|i can be solved in terms of the 
method so-called "The iterative numerical process" which is 
introduced in literatures (for example, see We have 

improved this method, and then fix the parameters a, b, c, d 
by fitting the well measured spectra of bottomonia. In our 
previous work (U, we explain the reason for the choice of the 
input for charmonia. However, for bottomonia, the situation 
is slightly different and the masses of T(1S), Xfco(lP), 
Xbi(lP), Xb2(lP) and T(2S) are chosen for the obtaining the 
values of a, b, c and d. Similar to the procedure we took in 
our previous work |8y, instead of directly fitting the masses, 
we construct a series of relations which should be fitted: 



' m[T(2S)] 


-m[T(lS)] = 


E [2 


3 Si]- 




m[T(2S)j 


-m[xbo(lP)] 


= E 


2 3 Si 


- E 


l 3 Po 


' m[T(2S)] 


-m[x6i(lP)] 


= E 


"2 3 Sr 


-E 


'l 3 Pi' 


, m[T(2S)] 


-m[ X62 (lP)] 


= E 


'2 3 Sx" 


-E 


'l 3 P 2 ' 



(5) 



where, E [n^ s+1 lj] represents the eigen-values of the radial 
equations (|4]i with various quantum numbers n r , j, 1, and s 
for the bottomonia. Because the parameters a, b, c and d are 
involved in the potential (0, E [n^ s+1 lj] must be functions of 
these parameters, m [meson] are the masses of the individual 
states which are shown in the following Tab. UJ 19]. 

Sequentially, the parameters a, b, c and d are obtained by 
solving Eqs.((5]l. By employing the Newton's iterative method 
(Th e details about the numerical method can be found in Ref. 
COll .'). we have achieved: 

a = 1.2165, b = 1.2988, c = 0.8686, d = 0.5886 (6) 

Here we set a s — 0.284 and k = 0.42 GeV which seem 
somehow different from the values given in literature 1 15U2U 
l22ll . But as noticed, the deviation may be included in the phe- 
nomenological parameters a, b, c and d. The choice of a s has 
another reason which is associated with our treatment of the 
contribution of the Lamb shift (see next section), anyhow this 
value is not very far apart from that given in literature. 

Given a, b, c and d in (|6), the masses of the bottomonia 
states are determined as: 

M(n? s+1 l j )=£;[n I 2s + 1 l j ] + E (7) 
where, Eq is the zero-point energy: 

£? = m[T(lS)]-£[l 3 Si] (8) 
and the final results are shown in the Tab. [Jbelow. 



Explicitly, in the process, the masses of the mesons with 
superscript "fit" are taken as inputs to obtain the parameters 
and then the masses of other states in the table are predicted. 

For readers' convenience and a clear comparison, we also 
list the results for the charmonia which were obtained in our 
previous work [8J . 



TABLE I: The mass spectra for the bottomonia states (in GeV), 
with irib = 4.8GeV. The A/exp is the value of the mass given 
inPDGHl. 



meson 


Mexp 


Prediction 


meson 


Mexp 


Prediction 


^(PSo) 


9.3020 


9.4124 


r? 6 (2 1 S ) 




9.9932 


T(l 3 Si) flt 


9.4603 


9.4603 


T(2 3 Si) flt 


10.0233 


10.0233 


XM)(l 3 Po) fit 


9.8594 


9.8594 


T(l 3 Di) 


10.1611 


10.1607 


XH(l 3 Pl) flt 


9.8928 


9.8928 


T(1 3 D 2 ) 




10.1719 


MPPi) 




9.8897 


T(1 3 D 3 ) 




10.1827 


X i >2(l 3 P 2 ) flt 


9.9122 


9.9122 


T(3 3 Si) 


10.3552 


10.3949 



TABLE II: The mass spectra for the charmonia states (in GeV), with 
m c = 1.84GeV. The A/exp is the value of the mass given in 
PDGfl9|l. 

meson A/exp Prediction meson Mexp Prediction 
^(PSo) 2.9803 3.0189 Xc2(l 3 P2) Kt 3.5562 3.5564 
J/^(l 3 Si) flt 3.0969 3.0969 ?7 c (2 1 So) flt 3.6370 3.6370 
Xco(l 3 Po) flt 3.4148 3.4148 i/)(2 3 Si) 3.6861 3.6861 
Xci(l 3 Pi) flt 3.5107 3.5107 i/>(3 3 Si) 4.1164 
ft e ( llp i) 3.5259 3.5100 



with 

a = 1.1715, b = 1.2250, c = 0.8087, d = 0.5291 (9) 

IV. THE MASS SPECTRA OF THE BOTTOMONIA AS 
THE LAMB SHIFT IS TAKEN INTO ACCOUNT 

As well known, the Lamb shift is due to the vacuum fluc- 
tuation and may cause sizable effects on the meson spectra. 
Indeed, the QED Lamb shift may not be very significant be- 
cause of smallness of the fine structure constant a l23ll . in 
fact, it just reaches order of 10~ 7 eV, but for the QCD case, 
the situation is different. 

In the previous section the effects of the Lamb shift are not 
included in the eigen-energy ©. Thus in this section, we take 
the Lamb shift into account. However, we do not introduce 
the Hamiltonian induced by the Lamb shift into the differen- 
tial equation because the corresponding pieces are very com- 
plicated. Instead, according to the traditional method, we cal- 
culate the effects in terms of the wavefunctions obtained with 
the original Hamiltonian, i.e. accounting < 4 , |i?Lamb|^' >> 
where i^Lamb is obtained via the loop diagrams and simply 
add the estimated values into the binding energies of various 
states. Including the Lamb shift effects in the expressions of 
the spectra, we re-fit the data to obtain the parameters a, b, 
c and d again and predict the mass spectra of the rest reso- 
nances. 

Namely, we set the masses of a bound states to be the mea- 
sured values: 

2m b + E + AE LM = M EX p, (10) 

where E is the solution of the eigen-equation, A£?lm is the 
energy caused by the Lamb shift and Mexp has been defined 
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in Tab. HI and HI1 already. Solving the equation, one can obtain 
the parameters again. 

The authors of Ref. lfl2l - [l4tl gave the theoretical expressions 
for the binding energies which involve contributions of the 
Lamb shift. It is well known that the induced Hamiltonian 
contributing to the Lamb shift is due to the vacuum fluctuation 
and can be obtained by calculating the loop diagrams order by 
order. Thus the Lamb Shift starts at 0(a 2 )|[l3|] 2 . The Lamb 
shift is: 

AE LM =< *|VLamb|* >, 

where ^ is the solution of the Schrodinger equation contain- 
ing only the Coulomb piece, can be written as : 

A£ L M[n,j,l,s] = m\AE(a 3 s ) + AE(a 4 s ) + AE(a 5 s ) 



+AE(a 6 s ) 



(11a) 



For illustrating the contribution of the Lamb shift to the 
spectra, let us directly copy Titard's formulas 1 1311 below, 
where we dropped the tree-level terms and the relativistic cor- 
rections, and we have: 



AE(a 3 a ) 
AE(at) 



C 2 

-a£7r^(2A>7£+4ai): 



'87m 2 

r 2 

4 ^ F 



A), 2 



(ai +7£y 

A 2 



s 4n 2 7r 2 

8 ; v 12 1E> 4 



(lib) 

jE(ai/3 

(11c) 



where, the n in Eq.dTTb stands for the principal quantum num- 
ber as n = n r + 1, where, n r and 1 are defined in Sec. [Til] 
All the constants as ai i _a2 j _bi, (3i (i = 1, 2, 3) are given in 
Ref.JH (also see HI |25fl28ll). Further more, Hoang et al. 
estimated the contribution of higher orders up to 0(a 5 s ) and 
0(af ) to the binding energies(See UH]). 

When we calculate the QCD Lamb shift effects for the 
hadron spectra, the potential not only contains the Coulomb 
piece, but also the confinement, so that the expression would 
be more complicated than that shown in Eq.(llb) and (11c). 
In fact, there cannot be analytical expressions for the expecta- 
tion values. 

The Lamb Shift A^lmIii, j, 1, s] depends on the coupling 



2 We indicate that the Lamb shift effects start from 0(a 2 ), and it seems 
that this allegation conflicts with Eq.jTJJ. As noted, in the expression 
Eq.{TT} which is given in Ref.| 13], A_Elm is not the potential derived 
from the loop, but the expectation value of the potential with the wave- 
functions which are solutions of the Schrodinger equation containing only 
the Coulomb piece. Since such solutions possess an exponential factor 
exp —a s fir, the expectation value of any function of r as f(r) should be 
proportional to a™ (n > 1), thus the AElm in Eq.(ll) start from erf. 
But indeed the potential pieces corresponding to the Lamb shifts directly 
come from the loop diagrams which start from O (a' 2 s ) . 



constant a s in Eq.dTTI) as[ 13]: 

2tt 



OL s {H ) = 



/? ln/VA 
1 



1 - 



gi ln(ln/i 2 /A 2 ) 
ffi ln M 2 /A 2 

/3 2 ln 2 (bV/A 2 ) 



/^ln 2 p?jk 2 
-/3?ln(Iii/i 2 /A 2 )-i8?+|8aA)]}. (12) 

Using the formulas given above, one can evaluate the Lamb 
shift of the charmonia states. The choice of the renormaliza- 
tion point /1 is suggested b y Pi neda et al., and "a natural value 
for this parameter " is lfl3l, l24ll : 



na B 



where, 



OB 

2\ 



"sO ) 



{l + (oi 



(13a) 



(13b) 



a 



2 ' 7T 



1e ai/3o + 



+ ( h7p — 

42 1E1 4 



IT 2 )' 



(13c) 



In the expression of the newly derived Hamiltonian there is 
a term In 2a fir jr (after a Fourier transformation from the mo- 
mentum space to the configuration space), which is UV diver- 
gent. To deal with the divergence, it is suggested to take an ef- 
fective method. For smaller range of r the Coulomb piece 1 jr 
obviously dominates, so that in < ^\HLamb\^ > one can use 
the wavefunction ^ which is the solution of the Schrodinger 
equation containing only the Coulomb potential, namely we 
can have an analytical solution for this asymptotic situation. 
Thus < *| \n2a(ir/r\^ >oc ln(na^/2). To make the UV 
divergence vanish, the suggested renormalization scheme is 
to set [i — 2 /nas- Indeed, in 11311 . other three alternative 
schemes were also suggested, here we just take this one and 
find the value of a s determined with this scheme is closer to 
that adopted in early literature for calculating the spectra of 
bottomonia. 

The value of the parameter A is chosen as 0.2 GeV for 
bottomonia[24), and at this point, 



= 0.284 



(14) 



which is the value of a s we used in Sec. [HI] It is noted that 
a s is different for different quantum number n: 

< =1 = 0.24 , a^= 2 = 0.284 , c^ =3 = 0.316. (15) 

We will use the n-related a s value for evaluating the spectra 
of the radially excited states of bottomonia. 

Simply adding the Lamb shift to the total binding energy 
is like that we change the zero-point energy for each state. 
We still select masses of T(1S), Xbo(lP), Xbi(lP), Xb2(lP), 
T(2S) as inputs, and solve the equation (0 again as we did in 
last section. But the value of a s in (0 is taken as that given in 
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Eq.(fT5b which depends on n. The new solutions of a, b, c and 
d are: 

a (LM) = L4256j 6 (LM) = 1 3553) 

c (lm) = o 8077; rf (LM) = 0.6849, (16) 

where the superscript LM refers to that all the corresponding 
parameters are obtained as the Lamb shift being taken into 
account. 

With these solutions, our predictions on the whole family 
spectra of bottomonia are presented in Tab. Hill 



TABLE III: The mass spectra with the Lamb Shift (in GeV), where, 
the LM stands for the contribution of the Lamb Shift, M is the pre- 
dicted mass when the parameters are set as in Eq.dl6b and M' stands 
for M' = M + AE LM - 



meson 


Ai?LM 


M 


M' 


Mexp 


^(LSo) 


-0.1064 


9.5274 


9.4210 


9.3020 


T(l 3 Si) flt 


-0.1114 


9.5717 


9.4603 


9.4603 


Xw(l 3 Po) fit 


-0.0618 


9.9212 


9.8594 


9.8594 


Xw(l 3 Pi) flt 


-0.0620 


9.9548 


9.8928 


9.8928 


-0.0621 


9.9507 


9.8887 




X 62 (l 3 P 2 ) fit 


-0.0622 


9.9744 


9.9122 


9.9122 


?? b (2 1 So) 


-0.0549 


10.0451 


9.9902 




T(2 3 Si) flt 


-0.0561 


10.0794 


10.0233 


10.0233 


T(l 3 Di) 


-0.0412 


10.2139 


10.1727 


10.1611 


T(1 3 D 2 ) 


-0.0412 


10.2272 


10.1860 




T(1 3 D 3 ) 


-0.0412 


10.2399 


10.1987 




T(3 3 Si) 


-0.0379 


10.4380 


10.4001 


10.3552 



V. THE SPECTRA OF THE foc(foc) MESONS 

In this section, we further study the spectra of the bc(bc) 
mesons. Except the ground state B c , the other states of the 
be mesons have not been well measured yet [ 19], we cannot 



directly fit the parameters from data as what we do for the 
charmonia and the bottomonia. It is noted that the parameters 
for charmonia and bottomonia are not drastically apart and 
since the bc(bc) family lies between charmonia and the bot- 
tomonia, we may interpolate those parameters for the bc(bc) 
family, namely average the values for charmonia and botomo- 
nia to be that for the bc(bc) mesons (See the Tab. IIVI >. 



TABLE IV: The parameters for foe (or fee) mesons 



meson 


a 


fo 


c 


d 




cc° 


1.1715 


1.2250 


0.8087 


0.5291 


0.36 


fofo* 


1.2165 


1.2988 


0.8686 


0.5886 


0.284 


foe ( or bc ) 


1.1940 


1.2619 


0.8387 


0.5589 


0.322 



"the parameters of the charmonia can be found in 1811 
*the parameters of the bottomonia are given in (6) 



Since the values of the aa s and ca s are more useful for the 
calculation as discussed before, we re-define: 

A = aa s and C = ca s . (17) 

Thus the parameters for the bc(bc) mesons are: 

/ A bc = (A b + A c )/2; b bc = {b b + b c )/2 

\ C bc = (C b + C c )/2; d bc = (d b + d c )/2. U5J 

The difference of the quark masses (m c = 1.8GeV and 
m b = 4.8GeV) make the Hamiltonian (Q]) to possess a more 
complicated form[9]: 

H = H Q + H 1 + ... (19a) 

2 

Hq = ?- + m b + m c + V{r) + S(r) (19b) 
Hi = H sd + H si (19c) 



- (V'(r) S'(r)) + + • S + V 2 V(r)S b ■ S c 

1 < ,y, W_ V . (r) ^ P.T)W.T) _l gt .^ (19d) 



m b m c \ r J \ r 2 3 



H si = U±- 2+ ±- 2 -J^\^V{r) + -^\ 2 -V>{r).I 2 +[ P \V{r)-rV\r)\ 
8 m b m c / Am b m c yr 

+ 2 (V(r) - rV'{r)) p 2 + i \ ~V'(r) + V"(r) - rV w (r)J | 



(19e) 



I 

So the Schrodinger equation we need is: where, 

2T¥(r) = (F + ffi)*(r) = (E + m b + m c )*(r) (19f) = m b m c (2Q) 

m b + m c 
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With this equation and the concerned parameters, we can 
predict the spectra of the members of the whole bc(bc) family 
shown in Tab.lVl 



TABLE V: The predict of the spectra of the B c meson 



Quantity 


our predict" 


our predict 6 


EFG[29, 30] 


KWLC[31] 


l 3 Si - 


VS 


0.0526 


0.0448 


0.0620 


0.0548 


2 1 S - 




0.5923 


0.5843 


0.5650 


0.5863 


2 3 Si - 




0.5723 


0.5744 


0.5430 


0.5795 


3 3 Si - 


2 3 S 1 


0.4021 


0.4067 


0.3540 


0.3652 


l 3 Po - 




0.3709 


0.3885 


0.3670 






l 3 Po 


0.0492 


0.0392 


0.0350 




l 3 Pi - 


l 3 Po 


0.0495 


0.0411 






1 3 P 2 - 


1'P1 


0.0290 


0.0320 


0.0280 




l 3 Di - 


2 3 Si 


0.1341 


0.1509 


0.1910 




1 3 Z> 2 - 


l 3 _Di 


0.0176 


0.0253 


0.0050 




1 3 D 3 - 


1 3 D 2 


0.0170 


0.0261 


0.0040 





"with the parameter in Tab. lIVI 

^Considering the effect of the Lamb shift and the parameters are taken as 

[ED 



VI. CONCLUSION AND DISCUSSION 

In this work, we study the role of scalar potential to the 
spectra of charmonia, bottomonia and the bc(bc) family. Our 
strategy is that the scalar and vector potentials have different 
fractions which manifest in their coefficients (in the text, they 
are a, b, c and d for the Coulomb and confinement pieces re- 
spectively). By fitting some members of charmonia and bot- 
tomonia which are more accurately measured, we fix them. 
Since except B c , the ground state of the bc(bc) family, other 
states have not been well measured yet, then we interpolate 
the parameters for charmonia and bottomonia to determine 
the concerned ones for the bc(bc) mesons. With those pa- 
rameters, we further predict the mass spectra of the rest reso- 
nances of charmonia, bottomonia and the whole bc(bc) family. 
It is shown that unlike the QED case where the fraction of the 
scalar potential is very small and negligible, for the quarkonia 
where QCD dominates, the fraction of scalar potential is of 
the same order of magnitude as the vector potential. This is 
consistent with the conclusion of Ref. 03211 and is not surpris- 
ing. As we indicated that for the vector-like coupling theories 
QED and QCD, the scalar potential can only appear at loop 
level or is induced by non-perturbative effect (QCD only). 
Thus it should be loop-suppressed. However, for QCD, the 
coupling is sizable and the higher order contributions and the 
non-perturbative effects somehow are significant, so that one 
can expect the fraction of the scalar potential is large. 

Moreover, the Lamb shift is induced by the vacuum fluctu- 
ation and only appears at loop level, indeed its leading con- 
tribution is at 0(a 2 s ). Therefore for the QED case, it is hard 
to observe the Lamb shift (observation of the Lamb shift is a 
great success for theory and experiment indeed), however, for 



QCD the effects are not ignorable. It is shown 0] that the 
NLO QCD effects may exceed the LO contributions at some 
processes. By taking into account of the Lamb shift, we re- 
fit the model parameters and find they are obviously distinct 
from that without considering the Lamb shift. 

The results help us to get better understanding of QCD, 
higher order effects and especially the non-perturbative ef- 
fects. Even though it is only half-quantitative, it is an insight 
to the whole picture. 

In this work, we adopt the renormalization scheme as fi = 
2/nas lfl3ll . which determines the effective coupling a s . It 
is worth emphasizing that a s depends on the principal quan- 
tum number n and this is different from that usually used in 
literature. But for the ground state of charmonia and bottomo- 
nia, the values of a s are quite close to that appearing in the 
literature. 

The predictions on the bc(bc) family will be tested at LHCb 
experiments where a great amount of the excite states of 
bc(bc) will be produced. Comparing with the data, we will 
learn more about the QCD and structures of the "final" meson 
family. 

Actually, in this work, we only consider the Cornell poten- 
tial which is supported by the area theorem and commonly 
adopted in phenomenological studies where the spectra and 
wavefunctions of heavy mesons are involved. Indeed, there 
are some other proposals. For example, the authors of Ref. 
ll33ll use the harmonic oscillator model to deal with the con- 
finement and further consider the effects of open charm loop 
for higher excited charmonia states which may induce energy 
shifts and change decay widths. Moreover, some authors in- 
troduce a phenomenological form for the spin-spin interaction 
ll34ll which may also result in energy level shift. In this work, 
we just restrict ourselves at the quark level QCD motivated 
potential whose form is given in Ref.Jld, but we may fur- 
ther our studies on the coupled-channel scenario in our com- 
ing work. 

The contribution of the scalar potential to the hadron spec- 
tra was noticed by some authors [38] long time ago, and 
its importance was confirmed. In our present work, we re- 
emphasize its role and discuss the origin in comparison with 
the QED case. In terms of the newly achieved data on charmo- 
nia and bottomonia, we analyze the hadron spectra and gain 
all the concerned parameters. We also investigate the signifi- 
cance of the Lamb shift phenomenologically. Then we go on 
discussing the spectra of bc(bc) family within the same frame- 
work, the results will be tested in the future experiments 
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Appendix : Check the legitimacy of the approximation 

[2n(E -V(r)) 

We investigate and elucidate the legitimacy of the approxi- 
mation adopted in the text: 



V 



2fx(E - V(r)) 



through a few examples. 

In fact, such problems have been thoroughly discussed 
in literature and even written in textbooks, for example in 
Ref. 13511 and the relativistic corrections to the Cornell poten- 
tial can be found in Ref . li36h . Here we just re-do the numerical 
computation to convince our readers and ourselves about the 
legitimacy of the approximation adopted in the text because it 
is very important for getting the spectra. 

The 0-th order Schrordinger equation is: 



and in the perturbative method, the relativistic correction is: 

1 r o„ , ,2 



AE a 



dr r 



[p 2 Rno(r)\ 



4m 3 

A£ n=1 = 0.008984; 
AE n = 2 = 0.001460; 
AE n = 3 = 0.000558. 



(26) 



On the other hand, as we use the approximation: p 4 



2fi(E -V(r)) 
solution of Eq. ( [221 



in Eq.(l22l. 3 and then find the numerical 



(27) 



Ef =1 = 


-0.1223; 


Ef= 2 = 


-0.0300; 


Ef= 3 = 


-0.0130. 



If we define: AEi — Eq — Ef, then we have: 



k +V(r) 



* = E ^ 



With the relativistic correction, it will be: 



r 9 4 

^1 + v(r) 

2fi { ' 4m3 



* = Ei^ 



(21) 



(22) 



Note: here we ignore other irrelevant correction terms such as 
the L-S coupling etc. because we only concern the p 4 term 
and the approximation. 

First, let us use the Coulomb potential as an example be- 
cause there exists an analytical solution. 

For the Coulomb potential: 



V(r) = -- 



(23) 



We have the exact solution (the eigen-energy and the wave 
function): 



AEf =1 = 0.007265; 
AEf= 2 = 0.001249; 
AEf= 3 = 0.0002288. 



(28) 



From the Eqs.(r26l) and ( 1281 ). we can find that the error is 
near 0.001 GeV. The relative error is: 

= (Eg + A£") - Ef 
(Eg + AE n ) + Ef 

x^ 1 = 0.68%; 

= 2 = 0.35%; (29) 
= 3 = 1.28%. 

With the Coulomb potential the Schrodinger equation pos- 
sesses an analytical solution, so it is easy to see the error. 
However, for the Cornell potential, there is no analytical so- 
lution available, so that we need to use the numerical solution 
for our analysis. 

The Cornell potential is: 



Eg 



Rio(r) 



where, 



a 

4n 2 

2Kh~ Kr 

3 



l -K 



m 
2m 



(2-Kr)e 



-K 



m 
~2~'' 



,Kr 



27 



e 3 



Kr 



(24a) 
(24b) 

(24c) 
(24d) 



K 



-ma 



If m = 1.84 GeV, and a = 0.5 (these numbers are just taken 
for an illustration, but not for real physics), then: 



V(r) = — 
r 



br 



(30) 



where, we set a = 0.5, b = 0.2 and consider the second radi- 
ally excited state with n = 2 for an example. The numerical 
solution of Eq.(l2"TT) is: 



EH= 2 = 0.856872 



and the numerical solution of Eq.(l22l is: 



Ef' 2 = 0.918242 



So we have: 

AE' n = 2 



Ef 



n=2 



E, 



n=2 



0.06137 



(31) 



(32) 



(33) 



77111— 1 

Eq 



-0.115 



Eg= 2 = -0.02875 
Eg= 3 = -0.012778 



(25) 

here, we may use Eo as well, but for a clear comparison we use E\ instead. 
The error is not great. 
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In the numerical solution of Eq.(l2Tb we have the wave func- 
tion R n i(r), and with the perturbative method, the contribu- 
tion of the relativistic correction is: 

= 0.0669.. (34) 
Finally, from Eqs.d33b and ( f34t . we have the relative error: 

(E%= 2 + AE n = 2 ) - Ef= 2 

Even though the relativistic error seems large, in fact the error 
is only at order of a few of MeV. The QCD Lamb shift gener- 
ally results in a few of tens of MeV, thus the error brought up 
by the approximation seems not too serious. 

Here one important point should be clarified. Directly cal- 
culating Eq. (f34b . one would have an unrealistically large re- 
sult. The reason is obvious that unlike the analytical solu- 
tion for the Coulomb potential, the behavior of the numeri- 
cal solution near the zero point (r — > 0) is not appropriate 
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